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ABSTRACT 

The use of Type la supernovae (SNe la) as cosmological standard candles 
is a key to solving the mystery of dark energy. Improving the calibration of 
SNe la increases their power as cosmological standard candles. We find tenta- 
tive evidence for a correlation between the late-time lightcurve slope and the 
peak luminosity of SNe la in the B band; brighter SNe la seem to have shal- 
lower lightcurve slopes between 100 and 150 days from maximum light. Using 
a Markov Chain Monte Carlo (MCMC) analysis in calibrating SNe la, we are 
able to simultaneously take into consideration the effect of dust extinction, the 
luminosity and lightcurve width correlation (parametrized by Amis), and the 
luminosity and late-time lightcurve slope correlation. For the available sample of 
11 SNe la with well-measured late-time lightcurves, we find that correcting for 
the correlation between luminosity and late-time lightcurve slope of the SNe la 
leads to an intrinsic dispersion of 0.12 mag in the Hubble diagram. Our results 
have significant implications for future supernova surveys aimed to illuminate the 
nature of dark energy. 

Subject headings: distance scale - methods: data analysis - supernovae: general 



1. Introduction 



One of the most important problems in cosmology today is to sol ve the mystery of 



dark energy, the unkno wn reason for the observed cosmic acceleration (jRiess et al. 1 11998 



Perlmutter et al. I Il999l ). The discovery of cosmic acceleration was made using Type la 



Supernovae (SNe la), as SNe la can be calibr ated to be good cosmological standard candles 
( )PhiUipslll993l ; iRiess. Press, fc Kirshnerlll995l ). Improving the calibration of SNe la increases 
their power as cosmological standard candles. 
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A SN la is a thermonuclear explosion that completely destroys a carbon/oxygen white 
dwarf as it approaches the Chandrasekher limit of 1.4 M Q o This is the reason SNe la can be 
calibrated to be good standard candles. The first challenge to overcome when using SNe la 
as cosmological standard candles is properly incorporating the intrinsic scatter in SN la 
peak luminosity. The usual c alibration of SNe la reduces the intrinsic scat ter in SN la peak 
luminosity to about 0.17 mag (|Phillipall993l ; lRiess. Press, fc Kirshnerlll995l ). The calibration 
techniques used so far are based on one observable parameter, the lightcurve width, which 
can be parametrized either as Am 15 (decline in magnitud es in the £?-ba nd for a SN la in the 
first 15 SN la restframe days after .B-ban d maximum, see iPhillipa (119931 )). or a stretch factor 
(which linearly scales the time axis, see iGoldhaber et al.l (120011 )). The lightcurve width is 
associated with the amount of 56 Ni produced in the SN la explosion, which in turn depends 
on when the carbon burnin g makes the transition from turbulent deflagration to a supersonic 
detonation Jwheeleri booih . 



There are theoretical reasons to expect the existence of additional physical parameters 
that characterize the lightcurve of a SN la. For example, iMilne. The, fc Leisingl (120011 ) 
modeled the late light curves of SNe la, and found that the nonlocal and time-dependent 
energy deposition due to the transport of Comptonized electrons can produce a 0.10 to 
0.18 mag correction to the late-time brightness of SNe la, depending on their position in the 
Ami 5 sequence. 

In this paper we show that late-time lightcurves of SNe la can be used to tighten the 
calibration of SNe la, and reduce their effective (post-calibration) intrinsic scatter in the 
Hubble diagram. We describe our method in Section 2, and present our results in Section 3. 
We summarize and conclude in Section 4. 



Method 



2.1. Defining the sample 



We use the ligh tcurve fits obtained using the Joint SuperStretch code for SN la lightcurves 
JWang et al. ll2006h R to derive the lightcurve parameters (epoch and apparent magnitudes 
at maximum light and Amis). This is necessary for consistency in making color and Ami5 
corrections. 



1 The exception recently found bv lHowell et ajj (|2006l ) is super-luminous and can be easily separated from 
the normal SNe la used for cosmology. 

2 The B and V lightcurves were fit using no secondary bump. 
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We select a sample of "normal" SNe la (from a total of 135 available) that satisfy the 
following requirements: 

(1) B max V max ^ 0.5, 

(2) 0.015 < z < .1; 

(3) B lightcurve data extending beyond 100 days from maximum light. 



The requireme nt of B max —V max < 0.5 selects the "Branch normal" SNe la (jVaughan et al 



19951 ; lBranchlll998l ). which is standard in studying SNe la as cosmological standard candles. 
The requirement of 0.015 < z < .1 selects SNe la which are nearby yet distant enough so that 
they are not affected by peculiar velocities due t o cosmic large scale structure th at modifies 
the Hubble diagram of nearby standard candles ( Wang. Spergel. fc Turner! 1998h. which has 
the la r gest impact on t he ne a rest SNe la (s e e, e.g ., Zehavi et al. ( 19981 ); Cooray fc Caldwell 
j2006h : lHui fc Greene! J2006h : IConlev et all toom . 



This yields a sample of 13 SNe la. Of these, we exclude two more SNe: SN1999cc and 
SN1993h. SN1999cc has only one data point in the B lightcurve between 40 and 100 days 
from maximum light, insufficient to determine a slope independent of the lightcurve model. 
For SN1993h, the late-time lightcurve slope appears to depend on the lightcurve model, 
due to the appearance of a small secondary peak in the lightcurve data at ~ 80 days from 
maximum light. For the final set of 11 SNe la, the B lightcurves are very close to a straight 
line at > 60 days from maximum light in the SN la restframe, and there are a minimum of 
two data points at t > 60 days from maximum light. 

Our final sample consists of 11 SNe la that have lightcurve data at > 100 days from 
maximum light in the SN la restframe, and late-time lightcurve slopes insensitive to the 
lightcurve model. Table 1 lists the SNe la in our sample. We include all the data that are 
required to allow others to reproduce our results. Note that Abu an d N H are not required to 
reproduce our results [see Eqs.Q-Q], but are given for reference. 



2.2. Fitting the Absolute Magnitudes 

To convert observed apparent magnitudes to absolute magnitudes we assume a flat 
universe, with Q m = 0.3, £7a = 0.7, and H Q = 70kms _1 Mpc -1 . The B absolute magnitudes 
are 

Mt = B max - (1) 

where B max denotes the apparent magnitude at peak brightness in the B band, and the 
distance modulus \i is 

fx(z) =m-M = 51og[d L (2)/10pc]. (2) 
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Table 1: Nearby SNe la with well defined late-time lightcurves. 



SN name 


z 


B m ax 


AB max 


V m ax 


Ami 5 


bit 


A6 W 


Nit 


reference 


SN1990O 


0.031 


16.607 


0.035 


16.533 


0.961 


0.0160 


0.0008 


4 


H96b 


SN1990T 


0.040 


17.454 


0.086 


17.537 


1.298 


0.0152 


0.0028 


6 


H96b 


SN1990Y 


0.039 


17.674 


0.306 


17.494 


1.217 


0.0153 


0.0025 


7 


H96b 


SN1992bc 


0.020 


15.207 


0.009 


15.207 


0.960 


0.0147 


0.0002 


10 


H96b 


SN1993ag 


0.050 


18.298 


0.026 


18.081 


1.382 


0.0142 


0.0030 


3 


H96b 


SN1997dg 


0.033 


17.143 


0.033 


17.117 


1.281 


0.0177 




2 


J06 


SN1998ec 


0.020 


16.670 


0.053 


16.406 


1.074 


0.0164 


0.0033 


4 


J06 


SN1998V 


0.017 


15.912 


0.029 


15.738 


1.150 


0.0163 


0.0012 


5 


J06 


SN1999ef 


0.038 


17.524 


0.117 


17.418 


1.055 


0.0160 


0.0049 


3 


J06 


SN2000dk 


0.016 


15.637 


0.021 


15.543 


1.457 


0.0136 


0.0029 


4 


J06 


SN2000fa 


0.022 


16.026 


0.051 


16.055 


1.140 


0.0164 


0.0028 


4 


J06 



Note. — Note that H96b refers to Hamuv et al. ( 1996fJ ). and J06 refers to Jha et al. ( 2006 ). Nn refers to 
the number of data points in the lightcurve at t > 60 days from maximum light. 



For z < 0.1, the luminosity distance di,(z) is well approximated by 



dr(z) 



cz 



l + 5 d 



Qo)z 



(3) 



with q Q = Q m /2 — Q\. Since we are only considering very low redshift SNe la, our results 
are insensitive to the assumed cosmological model. 

Without any corrections, the dispersion of the SNe la in our sample in the Hubble 
diagram is 0.35 mag (see Fig.5(a)). This is due to extinction by dust, as well as the in- 
trinsic dispersion in t he luminosities of SNe la due to the physical diversity of these objects 
( iWoosley et al.l 120071 ). Our goal is to model these effects phenomenologically using observa- 
tional data, in order to calibrate SNe la into good standard candles. 



According to ICardelli. Clayton, fc Mathid (119891 1 . due to dust extinction, all true B 
magnitudes differ from the observed B magnitudes by 5B = B — B obs = Rg5(B — V), 
where S(B — V " ) = ( B — V) — {B — V) obs , and Rb is approximately constant. Following 



van den Berghl (119951 ). we use 1Z to denote the color correction coefficient parametrizing 
the combined effects of dust extinction and intrinsic SN col or variation. Assuming that 
B m ax — Vmax = for unreddened SNe la (IPhillips et al.lll999l ). the dust correction to M B l 



obs 
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is TIE, with E = B max - V m 



There is a well known correlation between SN la peak luminosity and Am 15 , th e decline 
in br ightness (in units of mas) 15 da ys (in the SN restframe) after peak brightness (jPhillips 
19931 ). Following IWang et al. I (120061 ). we parametrize the Amis dependence of M^ s with a 
piecewise linear function, £>(Ami 5 ) = a (Amis — 1.1) + al |Ami5 — 1-1 1. 



We use the Markov Ch ain Monte Carlo (MCMC) technique (see lNeill (119931 ) for a review), 
illustrated for example in iLewis fc Bridle! (120021 ). in the likelihood analysis. The MCMC 
method scales approximately linearly in computation time with the number of parameters. 
The method samples from the full posterior distribution of the parameters, and from these 
samples the marginalized posterior distributions of the parameters can be estimated. In the 
MCMC analysis, we fit M* to 



M f B u 



M + AM C 



M + TZE + a (Am 15 - 1.1) + al \Am 15 - 1.1| + Pb H , (4) 

where h\_ t (in units of mag/day) is the late-time lightcurve slope (in the SN restframe). In 
our MCMC code, we use 

2 



X 



N 

E 



ABr, 



(5) 



The MCMC algorithm allows us to explore the dependence of M^ s on all the parameters 
simultaneously. 

The Hubble diagram dispersion is a conventional measure of the intrinsic scatter of SNe 
la. It is given by the standard deviation of the apparent magnitudes at maximum light of 
the SNe la, B max , from the theoretical predictions given by 

m (z) = fx(z) + M, (6) 

where M is a constant that denotes the absolute magnitude of SNe la. It can only be deter- 
mined through direct measurements of the distances to the SNe la using Cepheid variable 
stars. Since M has no impact on the relative calibration of SNe la that we study in this 
paper, we allow M to vary in order to obtain the smallest scatter in the Hubble diagram, 
independent of the M values obta i ned in the MCMC likelihood analysis. This is the same 
approach adopted by I Wang et al. I (120061 ) using a different likelihood analysis technique. 



3. Results 



We have used the restframe fitted lightcurves to obtain slopes of the lightcurves at late- 
time, chosen to be between 100 and 150 days from maximum light. Since the lightcurves 
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are very close to straight lines at > 60 days from maximum light, the measured slopes are 
robust and not sensitive to the lightcurve model. Table 1 lists the late-time lightcurve slopes 
bu for the SNe la in our samplejj 

Our data consist of (B max , AB max , E, Am 15 , bu) for each of the SNe la in our sample. 
The parame t ers th at we vary in our MCMC likelyhood analysis are (M, TZ, a, af, (3). 



Wang et al. I (l200q i found that TZ = 2.59 ± 0.24, a = 0.78 ± 0.09, and af = -0.54 ± 0.19. To 
avoid obtaining unphysical values of a and af due to the small size of our sample, we impose 
a prior of al > —2 in our analysis. 

We obtain millions of MCMC samples, which are then appropriately thinned for sta- 
tistical independence. We do not present results on the "nuisance" parameter M; it is 
marginalized over in obtaining the probability distributions that we present. Fig.l shows 
the marginalized probability density functions (pdf) of TZ, a, af, and (3. Fig. 2 shows the 
corresponding 2-D joint confidence contours (68% and 95%) for these parameters. For com- 
parison, Fig. 3 and Fig.4 show the pdf and joint confidence contours for TZ, a, and af, with 
(3 = 0. Fig. 5 and Fig.6 show the pdf and joint confidence contours for a, af, and (3, for 
fixing 7Z = 4.1 (standard Milky- Way dust). In Fig.l, Fig.3, and Fig.5, the dotted lines are 
the likelihood functions. Agreement between the pdf and the likelihood function indicates 
that the MCMC chains have converged and hence provide reliable pdf's. Table 2 lists the 
best fit, mean, and standard deviations for TZ, a, af, (3. 

It is reassuring that we obtain the same value of TZ whether or not we include late-time 
lightcurve slope bu in our analysis. Note that the prior of al > —2 has little impact if we 
include bu in our analysis (see Figs. 1-2), but it has a very large impact if we only fit for 
(M, TZ, a, al). This indicates that the set of parameters (M, TZ, a, af) does not properly 
describe the data, given reasonable physical priors on af. 

Fig. 7 shows the Hubble diagrams for three cases, with the solid line indicating mo (z) = 
fi(z) + M in each case, where M is obtained by minimizing the rms of \B max (z) —mo(z)\ (the 
Hubble diagram dispersion). The three cases are: (a) No corrections, (b) B max corrected by 
subtracting AM cor = TZE + a (Am 15 - 1.1) + af |Am 15 - 1.1|, with TZ = 2.45, a = 2.01, 
and af = -2.00. (c) B max corrected by subtracting AM cor = TZE + a (Am 15 - 1.1) + 
af | Amis - 1-1 1 + Pbu, with TZ = 2.41, a = 1.89, af = -1.45, and (3 = 77.49 (mag/day)" 1 . 
The parameter values in (b) and (c) correspond to the best fit values from the MCMC 
analysis (listed in Table 2). We find that including the correlation between luminosity and 
late-time lightcurve slope in correcting M^ s reduces the Hubble diagram dispersion from 



3 The uncertainties in bu are estimated from fitting the data points in the lightcurve at t > 60 days from 
maximum light to a straight line. Note that Abu values are not used in our analysis, see Eqs.(T4])-([5]). 



2 2.2 2.4 2.6 2.8 



1.6 1.8 2 2.2 
a 




Fig. 2. — 2-D joint confidence level contours (68% and 95%) corresponding to Fig.l. 
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31 



Fig. 3. — Marginalized probability density functions (pdf) of 71, a, and at, for (5 — 0. 




2.2 2.4 2.6 2.8 2.2 2.4 2.6 2.8 

-1 .8 r — ' ^ 




1.9 2 2.1 



a 

Fig. 4. — 2-D joint confidence level contours (68% and 95%) corresponding to Fig.3. 
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Fig. 5. — Marginalized probability density functions (pdf) of a, and af, and f3, for fixing 
11 = 4.1 (standard Milky- Way dust). 




-2-10 1 

a' 

Fig. 6. — 2-D joint confidence level contours (68% and 95%) corresponding to Fig.5. 
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Table 2: Best fit, mean, and standard deviations for 1Z, a, at, (3. 



parameter 


best fit 


mean 


rms dev 


including 6 H : 


1Z 


2.407 


2.445 


0.141 


a 


1.890 


1.870 


0.109 


at 


-1.453 


-1.418 


0.244 


(3 


77.490 


80.925 


12.779 




not including bu 


(0 = 0): 




K 


2.454 


2.480 


0.123 


a 


2.007 


1.994 


0.051 


at 


-2.000 


-1.969 


0.030 


forcing 1Z = 4.1: 


a 


1.372 


1.371 


0.345 


at 


-0.583 


-0.621 


0.633 




112.364 


117.889 


43.939 



0.14 mag to 0.12 mag. For comparison, fixing 7Z — 4.1 (standard Milky- Way dust) gives 
a Hubble diagram dispersion of 0.22 mag. This indicates that the assumption of standard 
Milky- Way dust is not supported by the data. 

To study the effect of the measurement uncertainty in B max for our sample, we have 
repeated our MCMC analysis setting AB max to a constant. This gives similar results as 
above. We find that including the correlation between luminosity and late-time lightcurve 
slope in correcting M^ s reduces the Hubble diagram dispersion from 0.13 mag to 0.11 mag. 
This is similar to what we find when using the actual AB max from Table 1. 

Fig. 8 shows the absolute magnitude corrected for color and Ami5 only [case (b) of Fig. 7], 
versus bu, with its coefficient f3 given by case (c) of Fig. 7. This shows that the evidence for 
the correlation between luminosity and late-time lightcurve slope is tentative due to the small 
size of the sample. A much larger sample of SNe la with well measured late-time lightcurves 
is needed to firmly establish this correlation. 
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3.5 4 

l0 g( CZ CM B ) 

Fig. 7. — Hubble diagrams for three cases, with the solid line indicating m (z) — fj,(z) + M 
in each case, where M is obtained by minimizing the rms of \B max (z) — 7710(2)! (^ ne Hubble 
diagram dispersion), (a) No corrections, (b) B max corrected by subtracting AM cor = TIE + 
a (Am 15 - 1.1) + af |Am 15 - 1.1|, with K = 2.45, a = 2.01, and af = -2.00. (c) B max 
corrected by subtracting AM cor = TIE + a (Am i5 - 1.1) + af \Am 15 - 1.1| + (3b lt , with 
TZ = 2.41, a = 1.89, af = -1.45, and (3 = 77.49 (mag/day)" 1 . For comparison, fixing 
72. = 4.1 (standard Milky- Way dust) gives a Hubble diagram dispersion of a — 0.22 mag, 
much larger than the Hubble diagram dispersions in panels (b) and (c). 
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19.5 



~l 1 1 - 



n 1 r~ 



-i 1 r - 



-i 1 r~ 



-19 



18.5 



-18 



$ Mg-«»..= B„„-[«E+a(Am 16 -l.l)+a'|im, 6 -l.l|]-M 
M B =(-19.97±0.14)+/Sb ll , /3=77.49 

_. i i I i , i I i i i L_ 



012 0.014 0.016 0.018 

b lt (mag/day) 



0.02 



Fig. 8. — The absolute magnitude corrected for color and Ami 5 only [case (b) of Fig. 7], 
versus the late-time lightcurve slope bit, with its coefficient (5 given by case (c) of Fig. 7. 
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4. Discussion and Summary 



Improving the calibration of SNe la increases the power of SNe la as cosmological 
standard candles. We have shown the tentative evidence for a correlation between SN la 
peak luminosity and late-time lightcurve slope in the B band; brighter SNe la have shallower 
lightcurve slopes between 100 and 150 days from maximum light (see Fig. 8). Using this 
correlation, we are able to reduce the intrinsic scatter of SNe la (measured by the Hubble 
diagram residue) to 0.12 mag (see Fig. 7). 

We use the MCMC algorithm in our likelihood analysis; this allows us to explore all the 
parameters simultaneously and efficiently in calibrating the SN la peak luminosity for our 
sample of SNe la. We only used AB max in calculating % 2 , since including the errors in color, 
Am 15 , and late-time lightcurve slope introduces spurious pdfs in the MCMC analysis. Our 
results are not sensitive to the uncertainty in B max (setting AB max to a constant gives similar 
results to using the AB max from Table 1). The parameter set is (M, 71, a, af, (3), which 
contains the nuisance parameter and correction coefficients for color, Am^, and late-time 
lightcurve slope respectively. The pdfs we have found for (TZ, a, af, (3) are very close to 
being Gaussian (see Fig.l and Fig. 3), except for the pdfs for af (which are cutoff by the 
prior of af > —2). Late-time lightcurve slope bu appears uncorrelated with color (see Fig. 2), 
which supports bu being an intrinsic observable that correlates with the peak luminosity of 
SNe la. 

The apparent correlation between peak luminosity and late -time lightcurve s l ope i n 
the B band that we have found is consistent with the results by ICappellaro et al.l (119971 ). 
They studied the late-time lightcurves of five SNe la in the V band with the focus of mod- 
eling lightcurves at very late times (later than 150-200 days from maximum light), and 
found a similar trend in the peak luminosity and late-time lightcurve slope (which they 
parametrized with the d ifference in V magnitude from maximum to 300 days, see Table 1 of 
Cappellaro et~al1 Jl997h l 



The coefficients for correcting for the correlation between decline-rate an d peak luminos - 
ity that we have foun d are very different from previous work (se e Table 2 and iPhillipd (119931 ) ; 
Hamuv et all Jl996ah ; IPhillips et all Jl999h ; IWang et al.l tood ) ) . This is likely partly a con- 
sequence of the particular sample of 11 SNe la we have used for our analysis. This is not 
surprising, since previous work by vario us groups also found very different Amis corrections 
( IPhillips et al.lll999l ; IWang et al. 1120061 ). indicating that this correction is highly sensitive to 
sample selection. 



Note that the color correction coefficient we found was not sensitive to whether we 
correct for late-time lightcurve slope, and is consistent with the color correction coefficient 
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of TZ = 2.59 ± 0.24 found by iWang et al. I (120061 ). This is very different from a coefficient 
of TZ = 4.1 that is typical of Milky Way dust. We found that fixing TZ — 4.1 gives a 
Hubble diagram residue almost twice as large as the Hubble diagram residue found when 
allowing TZ to vary (see Fig. 7 and its caption). This strengthens the evidence that the color 
correction coefficient TZ of SN data cannot be explained by extinction by typical Milky- Way 
dust. The deviation of TZ from TZ = 4.1 could be d ue to the mixing of intrinsic SN la color 
variation with dust extinction (jConley et al] 120071 ). or variations in the properties of dust. 
The defini tive resolution of this issue will likely require NIR observa tions of a large number 
of SNe la farisciunas. Phillips, fc Suntzefj[20o3 : IPhillips et aDl2006h . 



Our results show the benefit of extending the observations of SN la lightcurves beyond 
100 days in the SN la restframe, with several data points between 50 to 150 days to enable 
a robust measurement of the late-time lightcurve slope. Our sample only contains 11 SNe 
la, limited by currently available data. The expansion in the scope of SN observations to 
obtain late-time lightcu rves is straightforward for survey s which cover the same areas in the 
sky in regular intervals (IWangil2000l ; IPhillips et al.ll2006l ). 



Our results have significant implications for the survey strategy of future SN surveys, 
such as those planned for the Advanced Liquid-mirror Probe for Astrophysics, Cosmology 
and Asteroids (ALPACA)Q; Pan-STARRS^ the Dark Energy Survey (DESjj; the Large Syn- 
optic Survey Telescope (LSST)0, and the Joint Dark Energy Mission (JDEM). 
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